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Abstract 

We discuss the computation of the third virial coefficient in polymer systems, fo- 
cusing on an additional contribution absent in the case of monoatomic fluids. We 
determine the interpenetration ratio and several quantities that involve the third virial 
coefficient for star polymers with 4 and 6 arms in the good-solvent regime, in the limit 
of a large degree of polymerization. 



1 Introduction 



In the dilute regime the osmotic pressure of a polymer solution can be predicted successfully 
by using the virial expansion, which we write as 



n=l 



where c is the polymer number density, p the weight concentration, M the molar mass of 
the polymer, T the absolute temperature, and and R the Boltzmann and the ideal-gas 
constant, respectively. The virial coefficients S„ depend on the degree of polymerization 
N and on the chemical details. However, in the good-solvent regime renormalization-group 
arguments [1-3] indicate that, for ^ oo, the ratios 

A-n+l = Bn+lRg (2) 

where Rg is the zero-density radius of gyration, approach universal constants A^^^^ that are 
independent of chemical details and depend only on the polymer large-scale structure. 

Much numerical and experimental work has been devoted to the calculation of the second 
virial coefficient B2. Results for the higher-order coefficients are instead rare, both experi- 
mentally and numerically. In recent years some numerical computations of the third osmotic 
virial coefficient for solutions of polymers of different architecture have been reported. [4-11] 
However, in essentially all works an incorrect expression for the third virial coefficient was 
used. The correct expression, which is valid for any fluid of flexible molecules, was derived in 
ref. [9] , and used to determine the universal ratio A3 for linear polymers in the good-solvent 
regime. Let us report here the result. Let us consider a molecular fluid in which each molecule 
is formed by units interacting by means of an intramolecular potential \4ntra(i"i! • • • , t^n), 
where ri, . . . ,rAr are the unit positions. Molecules i and j interact by means of an inter- 
molecular potential Vlnter that depends on the positions {va^} and {ra }■ Given a quantity 
O which depends on the coordinates of two molecules, we define a zero-density average {■)^ 
as 

{0)%^^,,, ^ i- 1 . . . dr«drf . . . drg) O eM-PVinUW'^) - /?Mntra({r(^)})], 

Q2^ J drS^ . . . rfr^i^drf . . . rfrSJ) exp[-pv,^,,,{{r^''>}) - (3V,^U{r^''>})]. (3) 

The meaning of this average is easily understood: we fix the position of the first unit of the 
two molecules to avoid irrelevant volume factors and average over all possible conformations, 
weighting each conformation with the intramolecular Hamiltonian only (which correspond 
to consider the zero-density limit). Analogously, given a quantity O that depends on the 
coordinates of three molecules, we define an average (2) (3): it corresponds to averag- 

ing over all possible conformations of the three molecules keeping the first unit of the three 
molecules fixed in r^^^ , ^ , rf ^ . In terms of these quantities we define 

I2 ^ / rfV2 (/l2)S,r,„ (4) 
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(5) 
(6) 

where fij = exp(— /JVinter) — 1 is the Mayer function. The third virial coefficient is then given 
by: 

= —h-T,. (7) 

This expression contains two terms: the ffist one, proportional to I3, corresponds to the 
usual term that gives the third virial coefficient in monoatomic fluids. In addition, there is 
a second term that is not present in monoatomic fluids and is related to the flexibility of 
the polymer molecule. This additional term was neglected in refs. [4-8,10], and thus the 
corresponding estimates of the third virial coefficient are incorrect. Bruns, [4] starting from 
a general expression given in Yamakawa's book, [12] derives the correct expression for S3, 
but then he neglects Ti in the numerical calculation, stating incorrectly that it can be shown 
that such term vanishes for hard-core systems. The derivation of B3 given in ref. [9] does 
not give a physical interpretation to the additional term Ti. Here we present a different 
derivation that follows the approach of ref. [5]. It clarifies the physical meaning of Ti and 
explains why this term is necessarily present and non- vanishing. 

For linear polymers, even though Ti does not vanish, its contribution is small. Indeed, 
the results of refs. [8,9] provide the estimates 

Al = 9.80 ±0.02, 

Al= lim ( --/3i?-M = 10.60 ±0.04, 

N~*oo \ 3 J 

lim TiR-^ = 0.80 ± 0.05. 

N-*oo 3 

Thus, the contribution Ti lowers the third virial coefficient only by 8%. 

In this paper, we extend the calculations of ref. [9] to regular star polymers in which 
/ branches of equal molecular weight are connected to a single branching unit. [13, 14] 
This type of polymers is particularly interesting. First, they have several technological 
applications. [13] Second, they show a quite different behavior depending on the number / 
of branches, interpolating between linear polymers and hard colloids. Here we shall focus 
on the cases / = 4 and / = 6 with the purpose of investigating the quantitative role of the 
additional term Ti for a different polymer conformation. For both values of / we find that 
the additional contribution Ti is small: our results show that, for large A^, Ti/B^ k, 0.065 
and 0.05 for / = 4, 6, respectively. Note that the relative importance of Ti decreases as / is 
increased, indicating that star polymers become increasingly more rigid as / goes to infinity. 

The paper is organized as follows. In Section [2] we give a new derivation of the expression 
([7]). In Section [3] we explain the model we use and the simulation method, while in Section H] 
we present our results and compare them with the existing literature. 
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h 
Ti 



j d^Vi2d?Vi^ (/l2/l3/23)o_ri2,ri3 ' 



,ri2,ri3 



(8) 
(9) 
(10) 



2 A New Derivation of the Third Virial Coefficient for 
Flexible Molecules 



Equation ([7j) was obtained in ref. [9] by first performing an activity expansion in the grand- 
canonical ensemble. A different derivation is reported in ref. [4]. None of these two deriva- 
tions gives any physical insight on the origin of the term Ti and indeed Bruns [4] concluded 
incorrectly that Ti = 0. Here we present a new derivation that clarifies that Ti vanishes only 
if the molecules are rigid, i.e. if the probability of each conformation is density independent. 

We consider the general case of molecules that have many different internal conformations, 
labelled by an index a. Each conformation has a Boltzmann weight p^, Va oc exp(— /^V^ntra), 
normalized so that = 1. For convenience, we assume that the number of conformations 

is finite, as it occurs in lattice models, but the results are clearly valid also in the general 
case in which there is an infinite number of conformations (it is enough to replace sums 
by integrals). Particles interact by means of a pairwise potential \^(r, a,/5) that depends 
on the relative distance r (we fix a reference point on each molecule) and on the internal 
conformations. As usual, we introduce the Mayer function 

/ij(r, ajOj) = exp[-/3\/(r, a^, a^)] - 1. (11) 

Let us now assume that the number fraction of each conformation a is fixed. Then, the 
fiuid can be seen as a multicomponent mixture of simple molecules. In this case the virial 
expansion can be written as [15] 

Z = 1 - - ^ X„a;/3J2,a/3 - y ^«^/32;7^3,a/37 + 0{(?), (12) 

where 

/2,a/3 = j rfVi2(r,a/3), (13) 

hap-r = y rf^ri2rf^ri3/i2(ri2,a/5)/i3(ri3,a7)/23(|ri2 - ri3|,/37), (14) 

and c is the number density. If the molecules are rigid the number fraction Xa is density 
independent and equal to the zero-density probability pa- Thus, if we define 

h = ^PaPphaP h = ^PaP(3Pfh,a(3-„ (15) 

a/3 a/37 



we obtain 



- - 2^'"y^^ + ^^'^^' ^^^^ 

which is the usual virial expansion, with Ti = 0. On the other hand, for fiexible molecules x^ 
is density dependent (if Xa were density- independent, single-molecule properties, for instance 
the radius of gyration, would not depend on density, which is clearly unphysical). To derive 
the c dependence of Xa we proceed as in ref. [16], Section II. We consider a quantity R^"^ that 
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assumes the value 1 if the conformation one is considering is the a one, and zero otherwise. 
Exphcitly, given a configuration /3, the value i?^"'* of on this configuration is 

|o if/3^a ■ ^^^^ 

By definition 

Xo. = (18) 

In order to compute the virial expansion of the right-hand side, we consider L molecules in 
a volume V and write 

^ ^ J2i3,,...i3^ I dn ■ ..dYLR^p^Pp, ■ --Ph ]\i<j[^ + - rj\](3i,(3j)] ^^^^ 

Expanding the numerator we obtain 

/3i ft, ft 

( ^ 2 r 5Z ^ft PftPftPft^2,/32ft + • • • 

Analogously, for the denominator we obtain 



(20) 



^^+{^^¥"^-^2 + ... (21) 
Then, in the thermodynamic limit, L,V oo sX fixed c = L/V, we obtain 

Xa=Pa + CPa E Pph^P " + O(c^). (22) 

In order to check the correctness of this result we have verified that the condition ^ a;^ = 1 
is satisfied by our expression. If we now substitute this result for Xa into Equation f|T2|) . we 
obtain for 

B3 = - ^PaPpP^h,aph,a^ + (23) 

0/37 

It is easy to recognize that 

Ti = j dhi2dhi3PaPf3p^fi2{ri2,aP)fi3iri3,a-f) - L 



2 

2 ) 



a/37 



= E^"^'3P7/2,a/3/2,a7 - /|, (24) 
0^97 
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so that we reobtain ([7j). Note that Ti can also be written as 



) 



2 



(25) 




This relation shows that Ti > and that Ti vanishes only if 



^^7-^2,^7 — P-yh,l3i 



(26) 
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for any pair a and (3. Therefore, Ti = only if YlyP-y-^'^^ot-r is independent of a. If this 
condition is satisfied, we have ^^^7/2,07 = h, so that the density correction that appears 
in Equation fl22l) vanishes. Therefore, Ti vanishes only if Xa does not depend on c, that is if 
the molecules are rigid. 

3 Model and Simulation Details 

We are interested in determining universal properties in the good-solvent regime, in the limit 
in which the degree of polymerization goes to infinity. We can thus use any model that 
captures the basic polymer properties. For computational convenience we consider the well- 
known self-avoiding walk (SAW) model on a cubic lattice. A star polymer with / branches is 
represented by / SAWs starting at a common point. A regular star is given by the center ro 
and by / branches rjj, i = 1, . . . , f , j = 1, . . . , Nf, such that |ro — rj^i| = 1, |rjj — rjj+i| = 1, 
and all lattice sites are visited at most once. The total number of monomers is A^ = fNf + 1. 
This model is well defined up to / = 6. For larger values of /, one should use a larger core 
or a model with soft interactions as in ref. [17]. 

We simulate the model by using different types of moves: 

(i) We consider pivot moves [18-20] applied to a single arm (analogous moves were used on 
the tetrahedral lattice in ref. [21]). These moves have been shown to be very efficient 
in simulations of linear polymers. [20] In star-polymer simulations they are not equally 
efficient since they are rarely accepted when the pivot is close to the center of the star. 

(ii) We consider cut-and-permute moves [22] applied to a single arm (see Figured]). They 
have been shown to be quite effective in simulations of polymers grafted to impenetrable 
surfaces and speed up the conformational changes close to the center of the star. 

(iii) We use local moves that involve moving one or two monomers of the walk. 

For / < 5 one can generalize the arguments given in ref. [20] to show that this algorithm 
is ergodic. No such proof is available for / = 6, though we expect that the combination of 
local and non-local moves makes the algorithm ergodic in this case, too. 

The virial coefficients are determined by using the hit-or-miss algorithm discussed in 



refs. [9,23]. 
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Figure 1: The cut-and-permute move applied to a star polymer with Nf = 8 and / = 3. 
One first chooses an arm and a pivot point P. The new arm is obtained by connecting the 
subwalk PE to the center C and then the subwalk CP to the translated point E. 

4 Results and Discussion 

We perform simulations of star polymers for / = 4 and f = Q with 50 < Nf < 2000 and 
50 < Nf < 1000 in the two cases, respectively. Since = fNf + 1, the total number 
of monomers is quite large; this should allow us to probe the universal large-A^ regime. 
Results for the constants A2, A3, and for A^ = —^I^R'^ are reported in Table [TJ In all 
cases the additional term Ti gives a small negative contribution. Quantitatively we find 
T1/B3 ^ 0.065, 0.05 for f = 4,6. This is consistent with the idea that star polymers are less 
and less flexible as / increases, so that we expect Ti/i?3 ^ as / ^ 00. 

The data reported in Table [T] show a systematic dependence on Nf, see Figure El and 
thus a proper extrapolation is needed in order to obtain the large- iV constants A2 and A^. 
We use the same procedure illustrated in ref. [8], fitting the data to 

A{Nf) =A* + aNj^ + hNj^\ (27) 

where A is a universal exponent whose best estimate is [24] A = 0.515 ±0.017 (other results 
are reported in ref. [25]). The exponent A2 is an effective one that takes into account several 
correction-to-scaling terms: as in ref. [8], we take A2 = 1.0 ± 0.1. Of course, the previous 
expression is only the leading part of an expansion in inverse (non-integer) powers of Nf. 
To monitor the role of the neglected terms, we have repeated the fit several times, each time 
including only the data satisfying Nf > Nf^-^i^. Stable results for A2 are obtained by using 
all data, while good fits of A3 require A^/,min = 100. The corresponding results are: 

_ / 9.979 ±0.009 / = 4, , . 

^2 - \ 14.174 ±0.016 / = 6; 
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Table 1: Estimates of A2, A^, and of A3 = -i/a^-*^. 







/ = 4 






/ = 6 




Nf 


A2 


^3 


A, 


A2 


^3 


As 


50 


10.899(2) 


52.48(4) 


49.75(3) 


15.280(2) 


113.43(6) 


108.79(6) 


100 


10.620(2) 


49.16(4) 


46.40(3) 


14.953(2) 


107.55(6) 


102.95(7) 


150 


10.497(2) 


47.64(4) 


44.98(5) 


14.809(2) 


105.10(6) 


100.33(6) 


250 


10.379(2) 


46.33(4) 


43.57(3) 


14.666(3) 


102.60(6) 


97.78(6) 


500 


10.257(2) 


44.98(4) 


42.24(4) 


14.515(3) 


99.99(6) 


95.45(6) 


1000 


10.173(2) 


44.13(4) 


41.48(5) 


14.416(3) 


98.38(8) 


93.72(8) 


2000 


10.115(2) 


43.47(4) 


40.82(4) 
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Figure 2: Third virial invariant ratio A^ for / = 4 and f — 6 versus N A — 0.515. We 
also report the extrapolation function ^43 + aNj^ , A = 0.515, determined in the fit. 
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r 39.56 ±0.16 / = 4, 

3 - \ 90.1 ±0.4 / = 6. ^ 

The errors we quote include the statistical uncertainty and the systematic error due to 
the uncertainty on A and A2. To compare with the literature it is useful to define the 
interpenetration ratio * = 2(47r)"^/^A2 and g = B^/Bl = A^/A^. We obtain 

^ ^ r 0.4480 ±0.0004 / = 4, 

\ 0.6364 ± 0.0007 / = 6; ^ ' 

, ^ r 0.397 ±0.002 / = 4, 
^ \ 0.449 ±0.002 / = 6. ^ ^ 

For comparison we quote ^* and g'* for / = 1: [9] 

= 0.24693 ±0.00013, (32) 
/ = 0.3240 ± 0.0007. (33) 

Other estimates of ^* and g* for / = 1 are quoted in refs. [9,25,26]. 

There are no numerical results for g*. Our estimates increase with / as expected, but, 
for f = 6, g* is still far from the hard-sphere value 5/8 = 0.625, which should be valid for 
/ — cxD (a discussion of the behavior of g* for large values of / is reported in ref. [27]). 
Our results for \E'* are in reasonable agreement with the numerical ones reported in the 
literature. For / = 4, refs. [28-30] quote = 0.46, 0.467, 0.453 ± 0.007, while for / = 6 
they quote = 0.64, 0.665, 0.63 ± 0.01. The results of ref. [29] are those that differ more 
significantly. Note, however, that in this work much smaller values of Nj arc used; moreover, 
no proper extrapolation is performed. Field theory results differ instead quite significantly, 
predicting "if* = 0.517, 0.798 for / = 4, 6, respectively. [31] Recent experimental results 
for star polystyrene in benzene are reported in refs. [32,33]. They quote 0.43 ^ ^ 0.46 
for / = 4 and ^ ^ 0.60 for / = 6, in reasonable agreement with our results. They also 
estimate the factor g. The results show a strong dependence on the molecular weight: the 
two samples with highest molecular weight My, give g ~ 0.40, g ~ 0.43 for / = 4 and 
g ~ 0.39, g ~ 0.50 for / = 6, with g increasing with My^. These results are close to our 
estimates, even though the experimental results apparently prefer somewhat larger values. 
Note that similar discrepancies are observed for linear polymers, see the experimental results 
cited in refs. [32-34] and references therein. Older results for ^* are cited in ref. [35]. The 
experimental estimates of the interpenetration ratio for star polystyrene in toluene show a 
significant dependence on the molecular weight My,: for the largest values of M^, experiments 
give ^' ~ 0.46, 0.55 for / = 4 and ^' ~ 0.65 for / = 6. These results are reasonably close 
to our estimates. The experimental values of ^ for polybutadiene in cyclohexane quoted in 
ref. [35], 0.42 ^ ^ ^ 0.47, are also consistent. 

In conclusion, we have shown that Ti is small but not negligible: if Ti is neglected the 
error is of 6.5% and 5% for / = 4, 6, respectively. Moreover, our results allow us to compute 
the osmotic pressure in the dilute regime with good precision (we expect the error to be of 
order of a few percent below the overlap concentration, see ref. [9]). We find: 

ri + X + 0.397X2 + ... / = 4, . . 

\l + X + 0.449^2 + ... / = 6, ^ ' 

where X = B2C. 
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